Neel and Valence-Bond Crystal phases of the Two-Dimensional Heisenberg Model on 

the Checkerboard Lattice 
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I use an improved version of the two-step density matrix renormalization group method to study 
ground-state properties of the 2D Heisenberg model on the checkerboard lattice. In this version 
the Hamiltonian is projected on a tensor product of two-leg ladders instead of chains. This allows 
investigations of 2D isotropic models. I show that this method can describe both the magnetically 
disordered and ordered phases. The ground-state phases of the checkerboard model as J2 increases 
are: (i) Neel with Q — (n,n), (ii) a valence bond crystal (VBC) of plaquettes, (hi) Neel with 
Q = ( 7r /2 ! tt), and (iv) a VBC of crossed dimers. In agreement with previous results, I find that at 
the isotropic point J% = Ji, the ground state is made of weakly interacting plaquettes with a large 
gap A m 0.67Ji to triplet excitations. 



Frustration-induced magnetically disordered phases in 
two dimensions (2D) recently have attracted substantial 
interest [l[ . Frustrated magnets are known to display un- 
conventional ground states with, in some cases, a large 
set of low-lying degenerate singlet excitations that are 
still not well understood. Among models of frustrated 
systems, the Heisenberg model on the checkerboard lat- 
tice (HMCL) has recently been intensively studied by 
various techniques @, H, 0, S S • This model is seen 
as a first step in the investigation of the 3D pyrochlore 
model. The emerging picture is that at the isotropic 
point {J 1 — J2), the HMCL spontaneously breaks the 
lattice's translational symmetry. The ground state is a 
singlet made of a collection of weakly coupled plaquettes 
with a large gap, A w 0.7Jx, to triplet excitations. Away 
from the isotropic point, the situation is less clear. There 
is no single method which can capture the full phase di- 
agram. 

In this letter, I introduce an improved version of 
the two-step density-matrix renormalization group (TS- 
DMRG) [lil which, as I will show, is very convenient in 
the study of the HMCL and other 2D frustrated models. 
This new version is based on using the two-leg ladder, 
instead of chains, as the starting point to build the 2D 
lattice. The main insight in using the two-leg ladder to 
construct the 2D lattice comes from large N predictions 
that frustration often induces ground states in which 
the translationaly symmetry is broken. In the strong- 
coupling regime of the disordered phase of S = 1/2 sys- 
tems, the system is made of a collection of singlets or 
plaquettes. This strong coupling regime cannot be de- 
scribed starting from independent chains which are gap- 
less. Starting from a single chain, small transverse per- 
turbations can yield a gap within the TSDMRG. But this 
gap is often small and it is difficult to obtain reliable ex- 
trapolations. The two-leg ladder does not present this 
problem. It does already present a large gap A rj 0.5 
even in absence of frustration. Coupled ladders natu- 
rally evolve toward the 2D Neel state as the number of 
legs increases. Hence, in principle, disordered and or- 



dered phases could be described within a two-leg ladder 
version of the TSDMRG. This suggests that the two-leg 
ladder is a more natural starting point to describe ground 
state phases of 2D antiferromagnets than the single chain. 

Additional insights into this idea came from my com- 
parative studj^ of coupled chains with half-integer and 
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integer spins [lfJJ . In Ref. when starting from single 
chains, I found that although chains with 5=1 display 
the Haldane gap, A « 0.4, they converge much faster to 
the Neel state than those with S = 1/2. Furthermore, 
when a frustration induced disordered phase is present, 
it can be much more easily found in the case 5=1. 
Hence, following the equivalence between the two-leg lad- 
der and the Haldane spin chain, suggested by the Affleck- 
Kennedy-Lieb-Tasaki construction 12 1, it would be bet- 
ter to adopt the two-leg ladder as the building block for 
two-dimensional lattices. 

I will now illustrate this idea in the case of the HMCL. 
Following the usual notation, the HMCL is given by: 
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where < i,j > represents nearest-neighbor sites and 
stand for next-nearest neighbors on every other plaque- 
tte. J\ is set as the unit energy. 

The TSDMRG with ladders is similar to the method 
with chains. So I refer the reader to Ref.jl, 0| for a com- 
plete exposition of the algorithm. Here, I will discuss 
only briefly the main points of the algorithm. I start by 
dividing the 2D lattice into two-leg ladders; the Hamil- 
tonian ([1]) is written as, 



H = 



E 

ladders 



Hi a 



dder 



Hi* 



(2) 



where Pliadder is the Hamiltonian of a single two-leg lad- 
der, Hi n t contains the inter-ladder part. In the first step 
of the method, the usual DMRG method is applied to 



2 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

J, 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
J, 



FIG. 1: (a) Spin gap of the two- leg ladder for J2 = (cir- 
cles), 0.5 (squares), 1 (diamonds), 1.1 (triangles up), and 2 
(triangles down) . (b) Ground-state energies as function of the 
system size for a two- leg ladder for J2 = 1 (circles), J2 = 2 (di- 
amonds), and for the 2D lattice for J 2 = 1 (squares), J2 = 2 
(triangle up). 



generate a low energy Hamiltonian of an isolated ladder 
of N x sites keeping mi states. Then m,2 low-lying states 
of the superblock states, the corresponding energies, and 
all the local spin operators are kept. These energies rep- 
resent the renormalized low energy Hamiltonian of a sin- 
gle ladder. The Hamiltonian J2]) is then projected onto 
the tensor product basis of independent ladders, 



ladders 



(3) 



where ^ladder is an eigenfunction of iJ , ladder- This yield 
an effective Hamiltonian, 



ladders 



Odadder 



Hi* 



(4) 



The resulting effective coupled ladder problem which 
is ID is studied again by the DMRG method in the trans- 
verse direction. The TSDMRG is variational, as the orig- 
inal DMRG method, the subspace spanned by the wave- 
functions of the form 'J? is a subspace of the full Hilbert 
space of Hamiltonian ([I}. Its convergence depends on 
mi and 7712, the error is given by max(pi, P2), where pi 
and pi are the truncation errors in the first and second 
steps respectively. m 2 fixes the energy band-width 8E. 
The method is accurate only when the inter-ladder cou- 
plings are small with respect to 5E. In the present sim- 
ulations SE ~ 4. Since for the HMCL the inter-ladder 
and intra-ladder are of the same magnitude, in principle 
this approach would be plagued by the same deficiencies 
the block RG method. But if the starting point is cho- 
sen so that the essential physics is already contained at 
the level of the ladder, the effective strength of the inter- 
ladder couplings will be small even if the bare couplings 



FIG. 2: Short-range correlations Ci (circles), C T (squares), 
Cd (diamonds) for the two-leg ladder for uncrossed (a) and 
crossed (b) plaquettes as function of J2 



are not. This is particularly the case of models with frus- 
tration in which the competing interactions largely can- 
cel each other in the strong frustration regime, yielding 
weakly coupled sub-clusters. 

The ground state properties of an isolated ladder can 
readily be obtained. I keep up to toi = 96 and N x = 16 
and I target spin sectors from S z = to S z = ±4 and 
used open boundary conditions (OBC). The maximum 
error is pi = 1 x 1CP 4 . There is a gap A for all values of 
J2 investigated between and 2. The finite size behavior 
of gaps for some typical values of J2 are shown in Fig[T] 
The case J2 = reduces to the usual two-leg ladder which 
has been widely studied in the literature [l3(. For J 2 = 0, 
A w 0.5 . As J2 increases, A has a non-monotonous be- 
havior. This suggests a rich structure which is revealed 
more clearly by the analysis of the correlation functions. I 
computed the following short-range correlation functions: 
the bond strength along a leg Cj u _ e = (Si i iSi + i i i) ti ^ c for 
uncrossed (u) and crossed (c) plaquettes, the diagonal 
correlation C du c = (S^iSi+i^)^^, and bond strength 
along the rungs C r — (Si,iSi,2). Note that I have intro- 
duced a second index to the local spin. These correlations 
are shown in FigJ^l Four regions can be identified: (i) re- 
gion I (rung dimers): < J 2 < 0.6, Q u a < 0, d u » C ic , 
Cd„, c > 0, C r < 0, and \C r \ > IQJ; the dominant spin- 
spin correlations are along the rungs. The ground state 
properties of the ladder in this region are identical to 
those of the unfrustrated ladder (J 2 = 0). (ii) Region II 
(plaquettes I): 0.6 < J 2 < 1, C^ c < 0, |C/J > \C lc \, 
Cd U:C > 0, C r < 0, and \C r \ < |CjJ; the physics is 
dominated by that of the isotropic point. At this point, 
the ground state is a collection of weakly interacting un- 
crossed plaquettes. Both Cd c and C/ c vanish at J 2 = 1. 
In this region, the local spin configuration is the same 
on all the uncrossed plaquettes as shown in Figj3fb) . 
(iii) region III (plaquettes II): 1 < J 2 < 1.3, Q u < 0, 
C lc > 0, C du > 0, C da < 0, C r < 0, and \C r \ < \C lu \; 
in this region, the ground state is again dominated by 
uncrossed plaquettes. But now the local spin configura- 
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FIG. 3: The four phases of the two-leg ladder: rung dimers 
(a), plaquette I (b), plaquette II (c), and crossed dimers (d). 



tions on two consecutive uncrossed plaquettes are images 
of one another by reflection with respect to a plane pass- 
ing through the middle of the crossed plaquette between 
them. Region IV (crossed dimers): 1.3 < J2, Ci u < 0, 
C dc < 0, C la = C r = C du w 0, and \C r \ < |CJ; 
the ground state is dominated by the crossed dimers on 
crossed plaquettes as shown in FigJ3Jd). The sketch of 
the spin structure corresponding to each region is sum- 
marized in Fig[31 Since I applied OBC, for a given size, 
there are two possible ground states depending on the 
plaquette pattern: (a) ucu...ucu or (b) cue... cue. In re- 
gion I, the configurations (a) and (b) have nearly the 
same energy. This is consistent with the fact that the 
translational symmetry is not broken. But in Region II 
and III, (a) has the lowest energy, since it has a larger 
number of uncrossed plaquettes. By contrast, in Region 
IV where dimer order is dominant, it is (b) that has the 
lowest energy. 

The 2D systems are obtained by applying the DMRG 
on H e f f in the transverse direction. I studied systems of 
size N x x N y = 4 x 6, 8x10, 12 x 14 and 16 x 18. I kept 
up to to 2 = 96 and used OBC. Inter-ladder interactions 
will have very different effects depending on whether they 
correspond to a magnetic regime or a disordered regime. 
I will first consider their effects on region II, which in- 
cludes the isotropic point. Recently, there have been a 
number of studies which strongly suggest that the physics 
of the 2D systems is identical to that displayed by the 
two-leg ladder. In other words, the ground state is es- 
sentially made of weakly interacting plaquettes. If this 
is the case, it means that the inter-ladder interactions 
will not strongly modify the ground state wave function 
of decoupled ladders. Fig[T] shows that the ground state 
energy and A remain very close to that of an isolated 
plaquette. Thus in the vicinity of J2 = 1, inter-ladder 
interactions do not strongly renormalize the properties 
of an isolated ladder which themselves are close to those 
of an isolated plaquette. The extrapolated gap is found 
to be A = 0.67Ji which is in good agreement with the 
prediction from exact diagonalization [5J . The same con- 
clusion is seen in FigQ] for region IV where the crossed- 




FIG. 4: (a) 2D gaps as function of the system size J2 = 
(circles), 0.5 (squares), 1 (diamonds), 1.1 (triangles up), and 
2 (triangles down), (b) Correlation function along the legs as 
function of the distance for J2 = 1.25. 



dimer ground state found for the ladder is also the ground 
state of the 2D lattice. In both cases, the wave function 
made of the tensor product of the wave function of single 
two-leg ladders is a good variational wave function for 
the 2D system. In each case, the ground state energy 
of the 2D system remains very close to that of individ- 
ual plaquettes or crossed dimers. This can be explained 
as follows: when ladders are brought together to build 
the 2D lattice, the dominant local correlations are C; u 
in region II and Cd a in region IV; during this process, 
magnetic energy cannot efficiently be gained. For region 
II, this is because the two neighboring plaquettes of an 
uncrossed plaquette in the direction of the rungs involve 
frustrated bonds. Hence the system prefers the original 
configuration to avoid increasing its energy. For region 
IV, C r is very small. The system cannot increase it when 
the ladders are coupled, because the spins are already in- 
volved in strong diagonal dimers. There is, however, the 
possibility to gain magnetic energy by forming Neel order 
along the direction of the diagonal bonds ( J2 direction) 
as suggested in Ref.[7||. This is unlikely, however, because 
once such a phase is reached, I do not see how the sys- 
tem could go to crossed dimers at larger J2. The action 
of Ji which act as frustration in this regime decreases 
as J2 increases. Hence once this hypothetical Neel phase 
along the J2 bonds is reached, there is no obvious mech- 
anism that could destroy it as J2 increases to yield the 
crossed-dimer phase as suggested in Ref. Q • Such a Neel 
phase would be favored only when J2 3> J\. I made 
rough calculations with J2 = 4, 8 and I found that the 
system remains in the crossed-dimer phase. The situa- 
tion is apparently identical to the J\ — J2 chain where the 
independent chains regime is only reached in the infinite 
J2 limit. 

The situation is very different for regions I and III. In 
region I, the dominant local correlation is C r ; when the 
ladders are brought together, magnetic energy can be 
gained by an antiferromagnetic arrangement along the 
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FIG. 5: Ground-state phases of the 2D checkerboard model 
as function of J2. Note that the phase boundaries are rough 
estimates taken from the phases of the two-leg ladder. 



rungs. This enhances the local antiferromagnetic order 
which exists along the legs and ultimately leads to a Neel 
order with Q = (tt,tt). This is seen in the vanishing of 
the spin gap for J 2 = and J 2 = 0.5 shown in FigUJa). 
This is in agreement with results for J 2 = from quan- 
tum Monte Carlo (QMC) simulations Q and large S 
analysis I find that the TSDMRG ground state en- 
ergy —0.6011 at J 2 = 0. is not in ver y g ood agreement 
with the QMC result -0.6699 of Ref.[lJ]. Despite this 
discrepancy, the TSDMRG is nevertheless able to repro- 
duce the low-energy behavior of the ordered phase. This 
is not in fact surprising. In the Resonating valence bond 
picture, the Neel state and its low energy excitations can 
be written as linear combination of a tensor product of 
dimers. The TSDMRG variational solution of Hamilto- 
nian|T]) which is a linear combination of the wave func- 
tions ^ has exactly this form. A similar analysis also 
applies for region III. Ci u is dominant in region II. But 
as seen in Figj^l Ci u has a minimum at J 2 = 1 and then 
increases. It becomes very close to C r when J 2 enters 
region III. Hence magnetic energy can be gained again 
through the rungs. Since the structure along the legs 
is not modified from Figl^d), the resulting wave vector 
will be Q — (7r/2,7r). This is seen in FigU]in the behav- 
ior of the spin-spin correlation function Ci(i) along the 
legs. Ci(i) displays a period of 4. The correlations be- 
tween the rungs (not shown) oscillate with q y — tt. Fig[S] 
presents a sketch of the different ground state phases of 
the HMCL as function of J 2 . I note that in Ref.[7|], a very 
similar phase diagram was suggested; the only difference 
with the TSDMRG phase diagram is the wave vector of 
the Neel phase between the plaquette and crossed-dimer 
phases. 

In summary, I have shown that the TSDMRG method 
can reliably be used to study the disordered phases with 



short correlation lengths of isotropic 2D models. In these 
phases, the system is a collection of dimers or plaque- 
ttes. This makes the two-leg ladder a very good starting 
point for a variational calculation. I showed that the ba- 
sic physics of 2D systems could already be read through 
short-range correlations of the two- leg ladders. This vari- 
ational calculation is able to predict reasonably magnet- 
ically ordered phases as well. In this work, I did not dis- 
cuss the question of low-lying singlet excitations within 
the gap. Targeting them will lead to large truncation er- 
rors and the calculations will become impractical. These 
excitations are naturally truncated out when they are not 
needed to form a target state. Finally, The same method 
could be applied to the Sutherland-Shastry, J 1 — J 2 or 
the Kagome models in 2D. It could also be applied to the 
pyrochlore lattice, provided that the Hamiltonian could 
be written in some form involving ID subsystems with a 
large gap. 

This work started during a visit at the Weizmann In- 
stitute. The author thanks E. Altman for hospitality. 
This work was supported by the NSF Grant No. DMR- 
0426775. 
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